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ABSTRACT 

We report on the aperiodic X-ray timing and color behavior of the accreting millisecond X-ray 
pulsar (AMXP) IGR J17511-3057, using all the pointed observations obtained with the Rossi X-ray 
Timing Explorer Proportional Counter Array since the source's discovery on 2009 September 12. The 
source can be classified as an atoll source on the basis of the color and timing characteristics. It was 
in the hard state during the entire outburst. In the beginning and at the end of the outburst, the 
source exhibited what appear to be twin kHz quasi periodic oscillations (QPOs). The separation Av 
between the twin QPOs is ~ 120 Hz. Contrary to expectations for slow rotators, instead of being 
close to the 244.8 Hz spin frequency, it is close to half the spin frequency. However, identification of 
the QPOs is not certain as the source does not fit perfectly in the existing scheme of correlations of 
aperiodic variability frequencies seen in neutron star low mass X-ray binaries (NS LMXBs), nor can a 
single shift factor make it fit as has been reported for other AMXPs. These results indicate that IGR 
J17511-3057 is a unique source differing from other AMXPs and could play a key role in advancing 
our understanding of not only AMXPs, but also NS LMXBs in general. 

Subject headings: pulsars: general - stars: individual (IGR J17511-3057) - stars: neutron - X-rays: 
binaries 



1. INTRODUCTION 

Low-mass X-ray binaries (LMXBs) are neutron 
star (NS) or black hole systems with low-mass (M 
< IMq) companion st ars. Out of th e nearly 200 
LMXBs known so far (|Liu et all 120071) . 13 are ac- 
creting millisecond X-ray pulsars (AMXPs), i .e., they 
have shown coherent millisecond pulsations (|Patrunol 
2010). The neutron stars in LMXBs are believed to 
be s pun up by accretion to millisecond p eriods (see, 
e.g., fBhattac harva fc van den Heuvell I1991D . but why 
only some LMXBs appear as AMXPs is still an open 
question. These systems can be studied through the 
spectral (color-color diagram - CD) and timing (Fourier 
analysis) properties of their X-ray emission. Based on 
the paths traced on the CD and associated variability, 
NS LMXBs are classified eith er as Z or atoll source 
(|Hasinger fc van der Klisl [1989). Correlated with the 
position of the source in the CD, the Fourier power 
spectra of the X-ray flux variations exhibit different 
variability components. Apart from coherent pulsations, 
the power spectra also exhibit aperiodic phenomena: 
broad components (noise components) and narrow 
com ponents (quasi pe riodic oscillations; QPOs) (see, 
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e.g.. Ivan der K lis 2006, for a review). 

The coherent pulsations at the NS spin frequency are 
thought to be due to hot spots formed by magnetically 
chann eled accreted matter (see, e.g., iPringle fc Reed 
1972). The origin of the aperiodic phenomena is poorly 
understood. They are presumably mostly associated 
with inhomogeneities in the matter moving in Kcplc- 
rian orbits in the accretion disk, or in the boundary 
layer. The timescale for matter orbiting in the strong 
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gravity region very near to the compact object is of 
the order of milliseconds (the dynamical timescale 
t = (r 3 /GM) 1 / 2 ~ 0.1 ms at distance r = 10 km 
from a compact object of mass M = 1.4M ). This 
strong gravity region can therefore be pr obed by the 
QPOs with millisecond timescales (see, e.g.. Ivan der Klisl 
2006) . The first millisecond phenomena were found with 
Rossi X-ray Timing Explo rer (RXTE): in Sco X-l. twin 
QPOs in the kHz range (jvan der Klis et ahl 119961 ) and 
in 4U 1728-34 similar QPOs and also burst oscillations 
(oscillations near the spin freque ncy v s that occur during 
type I X-ray bursts; see, e.g., Strohm aver fc Bildstenl 
2006 for a review). In many LMXBs, Av, the difference 
between the twin kHz QPO frequencies, was at the 
frequency of burst oscillations Vburst ■ A beat mechanism 
and V burst — v s were suggested by I Strohmaver et all 
1996) . This led to the sonic-p oint beat-frequency model 
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Observations of var iable Av in Sco X-l 
(jvan der Klis et aLlll996[ ) and later in other sources were 
inconsistent with the sonic-point beat-frequency mode l 
and led to a modified version (lLamb fe Millerl l200l 
Also, the relativistic precession model ( Stella et al.l 
1999) was proposed, in which v s plays no direct role 
in the formation mechanism of QPOs. Observations 
of LMXBs like 4U1636-53 (jvan der Klis et all Il996a| ) 
which exhibited Av ~ Vburst/ 1 ^, clinched by the discovery 
of kHz QPOs in SAX J 1808.8-3654 with Av ~ v s /2 
by iWiinands et al.1 (|2003f) . led to the proposal of new 
models, involving resonances . The relativistic resonance 
model (|Kluzniak et al.l 120041 ) and spin-resonance model 
(jLamb fe Millerl I2003D were proposed which allowed Av 
= v s and /or Av = v s /2. F or historical accounts se e 
Ivan der Klisl (2006, 2008) and lMendez fc Bellonil (I2007T) . 
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A model which can explain all the observations is still 
awaited. 

IGR J17511-3057 was discovered on 2009 Septem- 
ber 12 during galactic bulge monitoring by INTEGRAL 
(Bald ovin et alj 120091 ). It showed X -ray pulsations at 
244.8 Hz during RXTE Proportional Counter Array 
(PCA) point ed observations which es tablished its nature 
as an AMP (jMarkwardt et al.l l2009h . It al so exhibited 
type I X-ray bursts and burst oscillations (jWatts et al.1 
2009a. lAltamirano et al.l 2010a). A minimum companion 
mass estimate (assuming the NS ma ss = 1AM P) and 
orbit al inclination = 90°) is 0.13M Q f M arkwardt et all 
2009) and the up per lim it to the distance is 6.9 kpc 
( Altamirano et aLll2010af) . In this paper, we discuss the 
aperiodic variability of the AMP IGR J17511-3057. 

2. OBSERVATIONS AND DATA ANALYSIS 

We analyzed all the 71 pointed observations with the 
RXTE PCA of IGR J17511-3057 taken between 2009 
September 12 and October 6. Each observation covers 
one to several satellite orbits and contains up to 15 ksec 
of useful data for a total of ~ 500 ksec. We exclude data 
in which the elevation (angle between the Earth's limb 
and the source) is less than 10° or the pointing offset 
(angle betwee n source position an d pointing of satellite) 
exceeds 0°.02 (jJahoda et al.ll2006D . Type I X-ray bursts 
have been excluded from our analysis (excluding all data 
between 100 s before the rise and 200 s after the rise of 
the burst). 

2.1. Color Analysis 

The intensity and colors are obtained from standard-2 
mode data which has 16 s time resolution and 129 energy 
channels. The hard and soft colors are defined as the ra- 
tio of the count rates in the 9.7-16.0 keV/6.0-9.7 keV and 
3.5-6.0 keV/2.0-3.5 keV bands, respectively. The inten- 
sity is the count rate in the 2.0-16.0 keV energy range. 
To estimate the count rates in these exact bands for each 
of the five proportional counter units (PCUs) of PCA, the 
count rates are linearly interpolated between the corre- 
sponding energy channels. Then, the background sub- 
traction is applied to each band using the latest standard 
background model Q. The source intensity and colors are 
obtained at 16 s time intervals separately for each PCU. 
These intensity and color values are divided by the corre- 
sponding average Crab values calculated using the same 
method as described above, but obtained for the day clos- 
est in time to each IGR J17511-3057 observation. These 
Crab normalized values are then averaged over all PCUs 
and over time to obtain one mean value per observation. 
This method of normalization with Crab (|Kuulkers et al.l 
I1994D is based on the assumption that Crab is constant in 
intensity and colors and is used to correct for differences 
between PCUs and gain drifts in time. 

2.2. Timing Analysis 

The Fourier timing analysis is done with the ~ 125 
/is time resolution (Nyquist frequency of 4096 Hz) event 

1 See http://heasarc.gsfc.nasa.gov/docs/xtc/pcancws.html for 
more details 
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Fig. 1. — Intensity, hard and soft colors as a function of time. 
Each point corresponds to one observation and different symbols 
identify different groups as indicated. The inverted triangle refers 
to the first observation in the early rise of the outburst. 

mode data in the 2-16 keV energy range. The light curve 
of each observation is divided into continuous segments 
of 128 s, which gives a lowest frequency of 1/128 s ~ 
8xl0 -3 Hz. The power spectrum of each segment is 
constructed using the Fast Fourier Transform. These 
power spectra are Leahy normalized and averaged to get 
one power spectrum per observation. No background 
subtraction or dead time correction is done prior to this 
calculation. The Poisson noise sp ectrum is est i mated 
based on the analytical function of I Zhang et al.l (|1995l ) 
and subtracted from this average power spectrum. This 
power spectrum is then expr essed in source fra ctional 
rms normalization (see, e.g., Ivan der Klisl [2006) using 
the average background rate in the 2-16 keV band for 
the same time intervals. 

To improve statistics, we combine our observations 
into seven groups. The observations in a group are 
chosen to be close in time and have similar colors (see 
Section 3.1). The pulsar spike at 244.8 Hz is removed 
from each average power spectrum by removing all the 
data points in the 244-246 Hz frequency bins at full 
frequency resolution. The power spectrum of each group 
is fitted using a multi-Lorentzian function which is a 
sum of several L orentzians in the , v m n.x representatio n 
as introduced by B elloni. Psaltis. &: van der Klisl (f2002). 
In this representation, the characteristic frequency 
is v max = vq-\/1 + 1/(4Q 2 ) where u is the centroid 
frequency of the Lorentzian, and Q is the quality factor 
given as Q = i/ /FWHM where FWHM is the full 
width at half maximum. The Lorentzians are named 
Li with characteristic frequency Vi, where the subscript 
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Fig. 2. — Color-color diagram: each point corresponds to one 
observation. The inverted triangle refers to the first observation in 
the early rise of the outburst. 

T is used to identify the type of component. The 
components are : Lh - break frequency, Lh - hump, Llf 
- Low Frequency QPO, LhHz - hecto Hz QPO, L( ow - 
(perhaps) a low frequency manifestation of lower kHz 
QPO, L t - lower kHz QPO and L u - upper kHz QPO 
(|Belloni. Psaltis. fc van der Kudl2002l ). 

3. RESULTS 

3.1. Colors and Intensity 

Figure Q] shows the evolution of intensity and colors 
with time. The seven groups are indicated in the fig- 
ure. The source was first detected on 2009 September 
12 and reached its peak intensity on 2009 September 14. 
As the intensity decreases, both the hard and soft col- 
ors decrease with time. A similar correlated behavior 
has been observed previously in t wo other AMX Ps, XTE 
J1807-294 and XT E J1751-305 (jLinares et al.l 2005 and 
Ivan Straaten et aLl 2005. respectively) out of the thirteen 
known so far. In other AMXPs and many non-pulsating 
transient atoll sources, the hard color remains constant 
while the soft color changes in corr elation with intensity 
(see, e.g. Jvan Straaten et al.l 120051) . The decrease in the 
colors has a steeper slope after MJD ~ 55103. The data 
in this steep fall are represented in the CD in Figure [5] 
by the points with a value < 1.12 in the hard and soft 
colors. The light curve shows bumps at MJD ~ 55090, 
55093, 55098, 55102 and 55107 which are approximately 
traced by the hard color but not so closely by the soft 
color. The soft color shows considerable scatter during 
the entire outburst which obscures any short-term corre- 
lation with intensity. Large scatter in colors is also seen 
in other AMX Ps, b ut only in the tails o f their outbursts 
(jLinares et al.l 2005; Ivan Straaten et all 2005). 

3.2. Aperiodic variability 

The multi-Lorentzian fits to the power spectra of the 
seven groups are shown in Figure |31 with the best-fit 
parameter values listed in Table. 1. We report a 
conservative measure of the single trial significance (a) 
of all parameters which is given by P/P_, where P 
is the rms normalized power and P_ is the negative 
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TABLE 1 



Note. — Best-fit parameters viz. Characteristic frequency 
i/mai, Quality factor Q and Fractional rms. Also given are single- 
trial sig nific ance a of the Lorentzian feature defined as P/P_, see 
Section 13.21 tentative c ompon ent identification (Ident.) (based on 
scenario 1, see Section 13.3.21 ) of the 7 groups. Errors are at 1 
standard deviation. 

error on P and is calculated using A% 2 = 1. For a 
good fit (reduced \ 2 ~ 0.8-1.2, see Table 1) to the 
power spectra, 3-5 Lorentzians are needed. In some 
groups, the coherences of some components have been 
fixed □ as they are insufficiently constrained by the 
data. The fixing of these parameters does not affect 
the values of the other parameters. Groups 3-6 require 
three Lorentzians and show power spectra similar to 
each other. Interesting components at high frequencies 
are seen in groups 1, 2 and 7, which therefore require 
additional Lorentzians in the fit. In each of groups 1 
and 2, two high frequency QPOs (~135 and ^260 Hz) 
are seen at similar characteristic frequencies. In group 
7, it seems that the same two components appear but 
moved to much lower frequencies (~70 and ~180 Hz); 
the centroid frequency difference Au between the two 
simultaneous peaks is similar to that seen in groups 1 

2 When the broad noise components give a negative coherence, 
Q is fixed to (which means it is a Lorentzian centered at zero 
frequency). For the other fixed parameters, the value of the fixed 
parameter is chosen to be the best-fit parameter value obtained 
when all the parameters are free. 
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and 2 between 100 and 150 Hz. 

The average power spectrum of just the five 
brightest observations in the peak of the outburst 
(which are the last 5 observations of group 1, hereinafter 
group la) is shown in Figure [3J The 44.7 Hz component 
seen in group 1, in group la is resolved into two 
components with characteristic frequencies of 30.9 and 
56.9 Hz. Similar behavior is exhibited by group 2 if we 
use six Lorentzians instead of five in the fit: the 36.6 Hz 
component is resolved into two components at 26.0 and 
47.0 Hz. The best fit parameter values are given in Table 
1. The centroid frequencies for group la are 30.9±1.2 
and 56.6±1.8 Hz, and for group 2a (the six Lorentzian 
fit of group 2) they are 25.7±1.2 and 46.6±1.0 Hz. So, 
these two components are close to being harmonics of 
each other in both cases. The reason why they are 
not exactly harmonics could be that with time the 
components moved slightly in frequency, varying in 
strength in different ways. Averaging the power spectra 
can then mask an ex act harmonic frequency relationship 
(jMendez et al.l 119981) . We note that although not all 
components in groups la and 2a have P/P_ > 3, they 
all appear in both, statistically independent, groups. 

3.3. Identification of the components 
3.3.1. Low-frequency Complex 

The average power spectra we observe at low fre- 
quencies resemble those seen in AMXP s and other 
atoll sources in the EIS (see, e.g., Ivan der Klisl 
2006). These sources are known to exhibit corre- 
lations among the frequencies of the powe r spec- 
tral components (fWiinands fe van der Klisl 1999, 
IPsaltis. Belloni. fe van d er Klis 1999. Ivan Straaten et al.l 
2005). We use these correlations as a tool to attempt 
to identify the components in our source. Figure H] 
shows the relation between Lb and Lh frequencies 
as obser ved in many atoll sources kno wn as the WK 
relation (jWiinands fe van der Kl is 1999), along with the 
data points of our source. The two lowest frequencies in 
each group follow this correlation. So, they appear to 
be Lf, and Lh, respectively. Similarly, the correlation of 
the coherence Q with characteristic frequency is shown 
in Figure [5] for Lh along with that in our source. The 
data for all but one groups fit the correlation. The 
feature at 8.3 Hz of group 5 deviates as it has a high 
coherence of 3.3. This feature could be a low frequency 
QPO (Llf) as has been observed before in AMXPs 
(Ivan Straaten et aT] I2005D and other atoll sources (see, 
e.g.. Ivan der Klisll2006l for a review). 

Figure [5] shows the existing the correlations between 
the characteristic frequencies of the different components 
in atoll sources versus the L u frequency. Assuming the 
highest frequency component in our source to be L u , the 
data of the seven groups are plotted. From the match 
of groups 1 and 2 to other sources, the lowest frequency 
components appear to be L\, and Lh, respectively. 
For the other groups, contrary to indications from 
Figure the lowest frequency components are instead 
closer to the Lh and Lf ow tracks. The fractional rms 
amplitude of the individual low frequency components 
in all groups are similar to the observed values of 



the same components i n oth er atoll sources in EIS 
(see Ivan Straaten etail 120051 for a comparison with 
AMXPs). 

It is known that the correlations in AMXPs are 
often shifted relative to tho se in other atoll so urces. 
Shift factors upto 1 .6 (see Ivan Straaten et all 2005, 
iLinares et al.l 2005, iLinares et al.l 2007! for details) 
have been observed so far. Except for groups 1 and 2 
(where no shift factors are required for Lb and Lh), we 
would require a shift factor 2.05 ( obtained using the 
same method as Ivan Straaten et aLl 2005) to make our 
two lowest frequency components fall on the existing 
correlations if we interpret them as Lb and Lh- As there 
is no single shift factor that would make all the data 
points of our source fit in the existing scheme of corre- 
lations, we do not apply any shifts in this paper, but 
consider other possibilities for this discrepancy between 
the correlations in Figures 2] and [5] We first discuss 
scenarios for the broad low frequency components: 

[A] . The two lowest frequency components are Lb and 
Lh respectively. Figure [7] shows the plot of the frequency 
of the two lowest-frequency and the highest-frequency 
component in each group , identified as Lb, Lh and L u , 
respectively, as a function of time. These components 
are seen in all the power spectra and their frequencies 
do not change much in time, so it appears that we 
are seeing the same three components in all the power 
spectra. Yet only the Lb and Lh of groups 1 and 2 fall 
on the Lb versus L u and Lh versus L u relation. This 
might be because the L u from the fits of groups 3-6 are 
at lower frequencies than usual as only a broad, blended, 
component can be fit due to poor statistics, which does 
not allow us to resolve other components present, if any. 
However, such a single broad high frequency component 
is common in similar sources in the EIS, and these do 
show Lb and Lh on the standard relation. For group 7, 
a narrow L u is seen so in this case that explanation for 
the discrepancy does not apply. Together, this makes a 
reliable identification of Lb and Lh difficult for groups 
3-7. 

[B] . The two lowest frequency components are Lb and 
Lh in groups 1 and 2, and in the rest of the groups 
they are Lh and Lg ow respectively, as suggested by the 
correlations in Figure |H1 This does not support our 
conclusion of observing the same components in all the 
seven groups inferred from Figure [7J 



3.3.2. High frequency QPOs 

The power spectra of groups 1 and 2 are unusual, in 
that the frequencies of Lb and Lh are appropriate to the 
EIS, but presence of high-Q twin high frequency QPOs 
is more like what we would expect in the lower left 
banana (LLB) branch (without the very low frequency 
noise, VLFN). The fractional rms amplitudes of these 
high frequency components in groups 1 and 2 are 
similar to the observed rms values of the twin kHz QPO 
components in other atoll sources in the LLB. The 
second highest frequency QPOs of groups 1, 2 and 7 
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Fi g. 3. — Power spectra of the seven groups and groups la and 2a (using six Lorentzians instead of five as in groups 1 and 2, see Section 
13.21 1 with multi-Lorentzian fit functions (see Table 1 for the best-fit parameters). Group numbers are indicated. The pulsar spike has been 
removed before rebinning in frequency. 

form a line well above the Li ow versus L u relation as 
seen in Figure [5] It is interesting to note that this line 
is also above the Li versus L u relation of other atoll 
sources. If we use the factor of 2.05 derived earlier in 
Section 13.3.11 to shift the high frequency components 
in v u and vg, they do fall on the existing versus v u 
correlation, but then the corresponding low frequency 
components do not fall on their respective co rrelations 
(for works done earlier see Ivan Straaten et ail 2005 and 



3 It has been suggested previously by van Straaten ct al. (2005) 
that Li and L( ow are different components 



ILinares et al.l 2005). Hence, as discussed in Section. 
13.3.11 we do not rely on the shift factors to help us 
identify the high frequency components. There are a 
number of possibilities for the identification of the high 
frequency components in groups 1, 2 and 7. We discuss 
below the possible scenarios: 

[1]. The high frequency QPOs seen in groups 1, 2 
and 7 are the twin kHz QPOs. The components at 44.7 
Hz in group 1 and at 36.6 Hz in group 2 are L( ow as 
they fall close to the Li ow versus L u correlation. Lg ow 
is expected to be seen in the EIS. However, harmonics 
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Fig. 5. — Relation of coherence Q with characteristic frequency 
for in IG R J 175 11-3057 com pared to other atoll sources (gray 
points; from van Straatcn et al. 2005). 

of Li ow have not been observed before, nor have twin 
kHz QPOs been seen in t he EIS (with one exception, 
4U 1728-34: iMiriiari et all 2003). 

[2]. The high frequency QPOs seen in groups 1, 2 and 
7 are the twin kHz QPOs. The components at 44.7 Hz 
in group 1 and at 36.6 Hz in group 2 are LhHz- kHz 
QPOs are often accompanied by LhHz components and 
harmoni cs of LhHz were observe d once previously in 4U 
1636-53 (jAltamirano et al.l l2~008) . However, the putative 
LhHz components have frequencies that are too low 
compared to all other sources, and as mentioned above, 
twin kHz QPOs are not expected in the EIS. 

[3]. The highest frequency QPOs are L u and the sec- 
ond highest frequency QPOs are LhHz- The components 
at 44.7 Hz and 36.6 Hz are L^ ow as in scenarion 1. The 
2 nd highest frequency QPO in group 7 could be either 
LhHz or Lg ow ; this is not clear from the correlation. 
The drawback of this interpretation is that LhHz and 
Li ow have never been observed together. LhHz is not 
expected in the EIS as observed from previous works: 
as seen in Figure [6l there are no LhHz components 
observed at low L u values characteristic for the EIS (see, 
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Fig. 6. — Characteristic frequencies of all the power spec- 
tral components plotted versus the characteristic frequency of L u . 
The gray points are no n-pulsating LM XBs. The data are from 
lAltaniirano et ai.l 12001) . and for Cir X-l. lBoutloukos et al.l (2006'). 
The seven groups of IGR J17511-3057 (black points) are indicated. 
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Fig. 7. — Frequencies of the two lowest frequency components 
(L;, and L^) and the highest frequency components (L u ) of the 
seven groups in IGR J17511-3057 as a function of time. 

e.g., Ivan dcr Klis 2006). 

[4]. The two highest frequency QPOs in groups 1 and 
2 are LhHz features. At centroid frequencies 138. lj^; \, 

respectively 

(139.2 ± 3.4, 258.1 ± 18.7 and 113.5 ± 9.7, 263.5 ± 14.0 
for groups la and 2a), they are all approximate har- 
monics of each other and as discussed in Section 13.21 
the features might have moved in time masking an 
exact harmonic relation. As L u is not present in this 
scenario, we cannot use the correlations in Figure [6] 
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Fig. 8. — Ratio of measurements of Av = u u - vg and spin fre- 
quency v s as a functi on of v s in IGR J 17511-3057 and other NS 
LMXBs (gray points; Altamirano ct al.| 2010b). The step function 
is Au/us = 1 for Us < 400 Hz and Avjv a = 0.5 for v s > 400 Hz. 
The curved line corresponds to Av = 300 Hz. The four AMXPs 
are marked with arrows. 

to identify the other components. The components at 
44.7 Hz and 36.6 Hz are Lt ow as in scenario 1. The 
drawbacks of this scenario are that the LhHz features 
are not expected to be observed in the EIS and the 
LhHz and Li ow have not been observed together. Also, 
the two highest frequencies of group 7 cannot be identi- 
fied in this scenario, as they are not harmonically related. 

It is clear that every scenario has its own short- 
comings, which makes a reliable identification difficult. 
If we try to identify the components just by comparing 
the power spectra with those of other sources, the 
identification does not comply with the state as deduced 
from the CD pattern. The source has some unusual high 
frequency components seen in groups 1, 2 and 7 and an 
unusual triplet of broad low frequency components in 
groups 3-6. The possibility that the highest frequency 
feature is not L u cannot be ruled out, especially for 
group 7. No single shift factor such as proposed for 
other AMXPs can restore a full match to the other atoll 
sources. 

4. DISCUSSION 

The behavior of IGR J17511-3057 in the CD and the 
power spectra at first sight appears similar to what has 
bee n observed in othe r atoll sources and AMXPs (see, 
e.g, Ivan der Klisl [20136!) . From the color diagrams and 
the shape of power spectra of groups 3-6, the source 
appears to be an atoll source in the EIS. The other 
AMXPs h ave also be en classified as atol l sources (see 
IWiinandsl et al. 200 6, IWatts et ail 20 09a. ILinares et all 
2008. iKaaret et all 2003. iReig et all 2000). However, 
closer study of the power spectral components indicates 
that this source is peculiar and does not fit well in the 
scheme defined by other sources. As discussed in Section 
I3.3[ all the possible scenarios for fitting our source in 
this scheme have their own shortcomings. Therefore, 
none of the components can be identified with certainty. 

The scenario that appears to require the least number 
of additional assumptions is scenario 1, where the only 



peculiarity is that twin kHz QPOs appear in what 
is otherwise an ordinar y EIS; this was rep orted once 
before, in 4U 1728-34 fMigliari et all 120031 ). We note 
that if the two high frequency components are twin 
kHz QPOs, the measured differences Av between the 
centroid frequencies of high frequency QPOs in groups 



1, 2 and 7 are 112] 



144.9±i|;| and 104.2t&jj Hz, 



respectively, and 119^5 6 and lSOj^j^g Hz in group la 
and group 2a, respectively. These values are inconsistent 
with being close to the spin frequency v s of 244.8 Hz as 
would have be en expected f or this so-called slow rotator 
{v s < 400 Hz; iMiller et ail 1998). Instead, they are all 
consistent with half the spin frequency, which otherwise 
has only been seen in fast rotators (v s > 400 Hz). This 
can be seen in Figure [5] which shows the plot of lAvjv s 
as a function of v s for AMXPs and other atoll sources 
(the spin frequency is inferred from burst oscillations for 
these systems). The step function shows the historical 
distinction between the slow and fast rotators. It has 
been sugge sted t hat Av and v s are ( nearly) independent 
(jYin et all 2007, iMendez fe Bellonif 2007). The curved 
line represents a constant Av of 300 Hz. To make the 
AMXPs SAX J1808.4-3658 and XTE J1807-294 fit this 
curve, t hey would have to be s hifted up by a factor 
of -1.5 (jMendez fe Bellonil 12001 . The points of IGR 
J17511-3057 for groups 1, 2 and 7 would require a factor 
~2.5 to fall on this curved line. However, groups 1 and 
2 do not require this same factor to fit in the scheme of 
correlations seen in Figure [6j but rather a factor 2.05. 
We applied no shifts to any data in this paper as there is 
no single factor. Note that all four AMXPs in Figured] 



are consistent with cither Av = v s or Av = 
shifts 



's/2 



without 



Our results favor mo dels like the s onic-p oint and 
spm-rcsonance model (lLamb fe Milled 120031) a n d the 
relativistic resonance model ([Kluzniak etaTI l200l 
which predict that either Av = v s or Av = v s /2. The 
relativistic precession model (jStella et al.lll999l ) predicts 
that Av should decrease whe n v u increase s as well as 

2006). 



decreases (see Figure 2.14. in Ivan der Klisl 2006). Our 
results do suggest a low Av, however the value is almost 
a factor of two lower than the value expected from the 
model at the observed v u of IGR J17511-3057. 

In conclusion, IGR J17511-3057 is indeed a very pe- 
culiar and interesting source. If scenarios 1 and 2 apply, 
the properties of the source can be summarized as fol- 
lows: a) It exhibits kHz QPOs while in the EIS, b) In 
spite of being a slow rotator, kHz QPO frequency sep- 
aration is Av ~ v s /2 and c) It requires different shift 
factors to fall on the frequency correlations of LMXBs at 
different times. Clearly this source could play a very im- 
portant role in testing the existing models for the origin 
of QPOs. More observati ons of this sou rce with RXTE 
or perhaps ASTROS AT (|Agrawall 120(51 ) . when it goes 
into an outburst again, are necessary to understand the 
nature of the different components. If we could observe 
the source in different spectral states, the components 
exhibited and their frequency evolution would help in 
establishing their nature. Observations of the high fre- 
quency QPOs and their evolution as a function of time 
and spectral state are key to their reliable identification. 



REFERENCES 



Agrawal P., 2002, In: 34th COSPAR Scientific Assembly, vol. 34 

of COSPAR, Plenary Meeting 
Altamirano D., van der Klis M., Mendez M., et al., 2008, ApJ, 

685, 436 

Altamirano, D., Watts, A., Linares, M., Markwardt, C. B., 

Strohmayer, T., & Patruno, A. 2010a, MNRAS, 409, 1136 
Altamirano, D., et al. 2010b, MNRAS, 401, 223 
Baldovin C, Kuulkers E., Ferrigno C, et al., 2009, The 

Astronomer's Telegram, 2196, 1 
Belloni T., Psaltis D., van der Klis M., 2002, ApJ, 572, 392 
Bhattacharya D., van den Heuvel E.P.J. , 1991, Phys. Rep., 203, 1 
Boutloukos S., van der Klis M., Altamirano D., et al., 2006, ApJ, 
653, 1435 

Hasinger C, van der Klis M., 1989, A&A, 225, 79 
Jahoda K., Markwardt C.B., Radeva Y., et al., 2006, ApJS, 163, 
401 

Kaaret P., in 't Zand J.J.M., Heise J., Tomsick J. A., 2003, ApJ, 
598, 481 

Kluzniak, W., Abramowicz, M. A., Kato, S., Lee, W. H., & 

Stergioulas, N. 2004, ApJ, 603, L89 
Kuulkers E., van der Klis M., Oosterbroek T., et al., 1994, A&A, 

289, 795 

Lamb F.K., Miller M.C., 2001, ApJ 554, 1210 

Lamb, F. K, & Miller, M. C. 2003, arXiv:astro-p h/0308179 

Linares M., van der Klis M., Altamirano D., Markwardt C.B., 

2005, ApJ, 634, 1250 
Linares M., van der Klis M., Wijnands R., May 2007, ApJ, 660, 

595 

Linares M., Wijnands R., van der Klis M., et al., 2008, ApJ, 677, 
515 

Liu Q.Z., van Paradijs J., van den Heuvel E.P.J. , 2007, A&A, 
469, 807 

Markwardt C.B., Altamirano D., Swank J.H., et al., 2009, The 

Astronomer's Telegram, 2197, 1 
Mendez M., van der Klis M., van Paradijs J., et al., 1998, ApJ, 

494, L65+ 

Mendez M., Belloni T., 2007, MNRAS, 381, 790 

Migliari S., van der Klis M., Fender R.P., 2003, MNRAS, 345, L35 



Miller M.C., Lam b F.K., Psaltis P. , 1998, ApJ, 508, 791 
Patruno, A. 2010, larXiv: 1007.11081 
Pringle J.E., Rees M.J., 1972, A&A, 21, 1 
Psaltis D., Belloni T., van der Klis M., 1999, ApJ, 520, 262 
Reig P., Mendez M., van der Klis M., Ford E.C., 2000, ApJ, 530, 
916 

Stella L., Vietri M., Morsink S.M., 1999, ApJ, 524, L63 
Strohmayer T., Giles B., Jahoda K, Lee U., 1996, In: Bulletin of 

the American Astronomical Society, vol. 28 of Bulletin of the 

American Astronomical Society, 1425 — h 
Strohmayer T., Bildsten L., 2006, New views of thermonuclear 

bursts, 113-156 

van der Klis, M., Swank, J. H., Zhang, W., Jahoda, K., Morgan, 
E. H., Lewin, W. H. C, Vaughan, B., & van Paradijs, J. 1996, 
ApJ, 469, LI 

van der Klis M., van Paradijs J., Lewin W.H.G., et al., 1996a, 

IAU Circ, 6428, 2 
van der Klis, M. 2006, Rapid X-ray Variability, In: Compact 

stellar X-ray sources, Cambridge Astrophysics Series, No. 39. 

Cambridge, UK: Cambridge University Press, 2006, p. 39 - 112 
van der Klis M., 2008, In: 37th COSPAR Scientific Assembly, 

vol. 37 of COSPAR, Plenary Meeting, 3283-+ 
van Straaten S., van der Klis M., Wijnands R., 2005, ApJ, 619, 

455 

Watts, A. L., Altamirano, D., Markwardt, C. B., & ., T. E. S. 

2009a, The Astronomer's Telegram, 2199, 1 
Watts, A. L., et al. 2009b, ApJ, 698, L174 
Wijnands R., van der Klis M., 1999, ApJ, 514, 939 
Wijnands, R., van der Klis, M., Homan, J., Chakrabarty, D., 

Markwardt, C. B., & Morgan, E. H. 2003, Nature, 424, 44 
Wijnands R., 2006, In: J. A. Lowry (ed.) Trends in Pulsar 

Research, 53 — h 

Yin H.X., Zhang CM., Zhao Y.H., et al., 2007, A&A, 471, 381 
Zhang W., Jahoda K, Swank J.H., Morgan E.H., Giles A.B., 
1995, ApJ, 449, 930 



